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PROBLEMS  STUDIED 

,  ‘>a 

■  -She  research  undertaken  under  the  present  grant  has  consisted  of  a 

number  of  investigations  into  nonlinear  wave  and  diffusion  equations.  These 
make  their  appearance  in  many  mathematical  modelling  problems,  and  we  have  been 
interested  in  applications  in  the  field  of  fluid  dynamics,  heat  and  mass 
transfer,  and  superconductivity.  The  problems  were  tackled  using  the  expertise 
available  both  through  permanent  staff  members  and  through  a  visitor 
programme,  and  equal  emphasis  was  laid  on  analytical  techniques,  numerical 

\  '  r  \ 

techniques,  and  modelling. 


SPECIFIC  RESULTS 
A 1 )  Lens  and  Antenna  design 

The  basic  problem  here  (originally  stimulated  by  a  query  from  the 
United  States  Navy)  is,  given  two  points,  to  design  an  optical  lens  which  will 
have  the  property  that  it  focuses  all  rays  of  light  from  the  one  point  onto 
the  other.  The  problem  reduces  to  the  study  of  an  unusual  pair  of  functional- 
differential  equations,  and  as  part  of  research  under  the  grant  it  was  shown 
that  if  the  lens  is  symmetric,  i.e.  if  the  two  faces  of  the  lens  are  the  same, 
then  there  exist  precisely  two  lenses  which  will  give  the  required  focusing, 
one  in  which  the  rays  of  light  do  not  cross  the  axis  of  the  lens,  and  the 
other  in  which  they  do. 

There  are  many  extensions  of  this  problem,  for  example  to  lenses  which 
are  not  symmetric  or  to  lenses  with  variable  refractive  index.  Investigation 
into  these  extensions  has  been  a  continuing  part  of  the  research  programme. 

One  general  aspect  of  which  has  been  studied  intensively  is  the  question  of 
how  well  a  slightly  imperfect  lens  can  perform.  In  several  situations,  there 
are  "best  possible"  lenses  and  these  include  off-axis  or  axially  displaced 
object  points,  and  lenses  focusing  at  two  different  wavelengths. 
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Or' ( 2 )  Blow-up  and  Quenching  for  Nonlinear  Diffusion  Equations 
u'  Equations  of  the  general  type 

ut  -  Au  +  f(u)  (1) 

can  exhibit  blow-up  solutions  (when,  at  some  finite  time  T,  u(x,T)  becomes 
infinite  for  some  value  of  x) ,  or  quenching  of  solutions  (when  at  some  finite 
time  T  a  positive  solution  becomes  identically  zero).  One  of  the  most  basic 
questions  to  be  answered  is  whether,  if  blow-up  does  occur,  it  occurs  at  just 
one  value  of  x  (or  perhaps  a  finite  number  of  values  of  x)  or  alternatively  at 
all  x  simultaneously.  As  part  of  the  research  under  the  grant,  it  was  shown 
that,  under  very  general  conditions,  at  least  in  one  dimension,  blow-up  does 
occur  at  most 'at  a  finite  number  of  of  points,  and  this  leads  to  a  number  of 
supplementary  questions  which  have  also  been  looked  at  under  the  programme. 

For  Example,  what  is  the  nature  of  the  profile  of  the  solution  just  before 
blow-up  or  quenching?  What  is  the  effect  of  introducing  spatial ly-dependent 
'coefficients  in  (1)?  Does  the  one-dimensional  result  extend  to  higher 
'  dimensions?  At  least  partial  answers  to  the  first  two  questions  have  been 
given.  The  third  question  remains  openj. 

3 )  Wedge  Entry  ^ 


This  is  the  problem  of  the  entry  of  a  wedge  into  water,  and  models  the 
impact  of  seas  on  a  ship's  hull.  The  investigation  arose  because  of  a  query  by 
the  UK  Navy.  Although  the  governing  equations  are  almost  classical,  little 
effective  work  had  been  done  on  them  beyond  numerical  calculations.  As  part 
of  the  research  programme  we  have  been  able  to  give  a  complete  account  of  the 
existence  and  qualitative  behaviour  of  the  solutions,  including  an  analysis  of 
the  asymptotic  behaviour  as  the  semi-angle  of  the  wedge  approaches  hm ,  i.e.  as 
the  wedge  deadrise  angle  becomes  very  small,  which  is  the  important  case  for 
naval  applications.  Additional  work  on  the  asymptotic  behaviour  of  slamming 
bodies  of  more  general  cross  section  has  also  been  carried  out  by  using  a 
"displacement  potential"  to  reduce  the  problem  to  a  variational  inequality. 
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(4)  Asymptotics  Beyond  All  Orders^. 

'Y  In  quite  a  number  of  recent  problems  situations  have  arisen  where 
asymptotic  methods  (say  as  a  parameter  c  -»  0)  give  the  wrong  answer  if  only 
algebraic  terms  in  c.  are  considered,  and  when  it  is  necessary  therefore  to 
take  account  of  terras  of  order  e'1/c.  Two  such  instances  are  an  equation 
arising  in  dendrite  growth, 
c  6  +  $  -  cos  9  . 

and  an  equation  arising  in  wave  propagation, 

.  £  v  +  v-l-v2. 

-'''A  systematic  investigation  of  such  problems  has  been  carried  out,  involving 
both  analytical  and  numerical  techniques. 

Higher-order  Diffusion  ^ 

i  Problems  involving  the  spread  of  viscous  droplets,  and  also  certain 
modelling  problems  in  semiconductors,  lead  to  higher-order  diffusion  equations 
of  the  type 


3fz  +  -o. 


.♦!«_  (  I  „  I  n  +  1  ’ 

3xz 


3x 


(2) 


The  droplet  case  is  the  case  m  -  1,  and  the  semiconductor  m 


We  have 


^carried  out  an  extensive  investigation  of  the  existence,  uniqueness  and 

^'qualitative  behaviour  of  similarity  solutions  of  (2)  in  the  parties' ar  case 

f 

'm  -  1  and  are  extending  the  work  to  other  values  of  m. 

V  S* 

tlacroscopic  Models  for  Solidification  ^  ~ 

The  mathematical  regularisation  of  ill-posed  Stefan  models  for 
supercooled  solidification  has  been  studied  from  several  aspects.  The 
introduction  of  surface  tension  poses  an  interesting  theoretical  challenge  and 
several  conjectures  have  been  made  concerning  the  resulting  morphology.  A 
regularisation  is  also  possible  using  a  "mathematical"  nucleation  to  create  a 
new  phase  boundary  in  such  a  way  chat  Che  we  1 1 - posedness  of  the  model  is  ensured. 
Yet  another  possibility  is  the  introduction  of  kinetic  undercooling  at  the 
free  boundary  and  this  powerful  mechanism  has  been  studied  analytically  and 
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asymptotically . 

Certain  explicit  solutions  to  two-dimensional  classical  Stefan  models 
have  been  catalogued  and  an  asymptotic  analysis  of  a  simple  continuous  casting 
machine  has  been  possible  in  the  typical  case  when  the  liquid  steel  from  the 
tundish  forms  a  long  thin  finger  in  the  caster. 

As  usual,  members  of  the  faculty  and  students  from  Oxford  have 
participated  in  tne  triennial  international  free  boundary  conferences,  most 
recently  held  in  Irsee,  Bavaria  (June  1987)  and  Montreal  (June  1990). 

''  ( 7 )  Mathematical  Models  for  Fibre  Tapering 

A  simple  hyperbolic  model  for  the  unsteady  axial  tapering  of 
cylindrical  optical  fibre  has  been  proposed  and  analysed  to  show  that  breakage 
is  only  possible  in  exceptional  circumstances.  This  work  has  spurred  further 
research  into  tapering  of  non-cylindrical  fibres,  and  in  particular  the 
problem  of  core  displacement  when  several  fibres  are  being  spliced  by  using 
^/"the  tapering  mechanism. 

^^V(8)^>Semiconductor  Fabrication 

The  derivation  and  study  of  hierarchies  of  nonlinear  diffusion 
equations  which  arise  in  modelling  semiconductor  dopant  diffusion  were 
initiated  during  the  early  part  of  the  contract  and  research  on  them  has 
continued  at  other  universities  ever  since.  The  most  recent  emphasis  has 
been  on  nonlocal  effects  resulting  in  pairing  of  vacancies  and  dopant  atoms: 
this  results  in  the  appearance  of  effective  diffusion  coefficients  which  depend 
on  the  boundary  data. 


><9)  Superconductivity  Modelling  .  ^ 

cy 

Work  has  been  initiated  on  the  analysis  of  macroscopic  models  for  phase 
transitions  in  superconductors.  The  Ginzburg-Landau  theory  (based  on  the  BCS 
theory)  leads  to  a  complicated  system  of  nonlinear  diffusion  equations  for  the 
magnetic  field  and  the  complex  order  parameter  which  describes  the  extent  to 
which  the  transition  has  been  completed.  It  appears  that  in  various 
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asymptotic  limits  this  system  can  be  reduced  either  to  vector  Stefan  problems 
or  to  transition  layer  problems  for  nonlinear  ordinary  differential  equations, 
depending  on  the  observer's  length  scale.  The  parameters  appearing  in  these 
reduced  models  characterise  boch  the  "type"  of  the  superconductor  and  the 

stability  of  the  intermediate  state. 

•*.  '  r 

» .( 10)  Hypersonic  Flow  ,  '  .  i  < 

’  ”  O'  J 

A  study  has  been  made  of  the  hyperbolic  system  describing  gas  flow  in 
the  inviscid  shock  layers  which  occur  in  the  vicinity  of  pistons  and  blunt 
bodies  when  Che  free  stream  Mach  number  is  large  and  the  specific  heat  ratio 
near  unity.  There  are  several  kinds  of  piston  motions  for  which  no  solution 
appears  to  exist  and  chese  situations  have  been  studied  using  a  regularisatior. 
technique  which  also  sheds  light  on  the  interesting  phenomenon  of  shock  layer 
/Separation  when  the  piston  is  abruptly  stopped  or  when  the  blunt  body  suffers 
an  abrupt  change  of  geometry.  . ''' 

--K11)  Combustion  in  Porous  Media  J 

Early  work  on  porous  medium  combustion  dominated  by  excess  gaseous 
fuel  has  been  modified  and  extended  to  combustion  in  the  presence  of  smaller 
amounts  of  gaseous  fuel.  The  type  of  stability  for  these  burning  tones 
terminated  by  temperature  or  combuscible  solid  exhaustion  has  been  determined. 
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